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Abstract 



We explore the equations of state and other properties of various quantum 
fluids (^He, ^He, their mixtures, and H2) confined within individual carbon 
nanotubes. Above a threshold number of particles, Na, the fluid density near 
the axis begins to grow above a negligibly small value. The properties of this 
axial fluid phase are sensitive to the tube size and hence to the transverse 
compression in the case of a bundle of nanotubes. We consider He at zero 
temperature and H2 at low temperatures. 



I. INTRODUCTION 

The physical properties of carbon nanotubes are expected to manifest a variety of in- 
triguing one- dimensional (Id) and quasi-ld behaviors We, and other groups, have 
found that fluids absorbed within these tubes ought to exhibit an "effective" dimensionality 
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which is sensitive to the geometry, species, number (A^) of particles and temperature (T). 
For example, in the case of isolated nanotubes, the ideal gas regime of very low A^, is pre- 
dicted to show crossover between Id and 2d behavior as T increases, due to the excitation 
of azimuthal motion of the atoms adsorbed on the inner wall of the tubes [|T^ . For a bundle 
of nanotubes, in contrast, it has been predicted and found experimentally that low gases 
which have small diameter are adsorbed strongly in interstitial channels, between tubes, in 
which the gases display Id motion |l5|J^, |TT| , H • Larger diameter gases adsorb preferentially 
within the tubes. [|l^,|l^] 

The novel collective properties of the absorbed interacting fluid at higher are now 
being explored. In this paper we consider the evolution, as a function of A^, of the structure 
and properties of dense quantum fluids within individual nanotubes at T = 0. Consider He 
for specificity. Above the ideal gas regime of density, there appears a regime of density in 
which liquid ^He is adsorbed on the walls of the tube; we call this a cylindrical "shell" phase 
r?! . At somewhat higher density this fluid should solidify; we expect that this film is similar 



to the much studied incommensurate monolayer solid film on graphite |[T8|-pO| . For values of 
N below a threshold value Na there exists only this shell phase. We predict here a transition 
in which there appears, for A^ > A"a, an axial phase of this fluid, signified by the presence of 
a fluid confined to the vicinity of the tube axis. This axial phase transition is qualitatively 
similar to two other transitions familiar in the field of adsorption: capillary condensation 



(CC) and the layering transition pTi|22| . While the latter (2d) transition is known to occur 
at finite T, the appearance of the axial phase is a genuine thermodynamic transition only at 
T = because the geometry is Id, so that fluctuations eliminate the transitions at higher T. 
This feature is also the case for a CC transition in a single pore, but often in experiments 
the effect of interactions with other pores is such as to permit a CC transition at nonzero 
T 0J^. Similar behavior of the axial phase transition should occur in a nanotube bundle 
("rope"), but we have not yet investigated that problem. Also, the same phenomenon is 
expected to occur for a classical system since nothing specifically related to quantum effects 
plays a paramount role |]^J2^ . 

We have investigated the axial phase with three distinct methods. The next section 
presents Hartree model calculations for four systems: "^He, ^He, their mixtures, and H2. 
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These calculations rely extensively on the use of other workers' calculations and data for 2d 
quantum systems pl|-p^ . Section III presents results of Path Integral Monte Carlo studies 
of H2. Section |V] reports results for ^He-'^He mixtures, obtained with the Hartree approach. 
Section |V| describes the results obtained with Density Functional calculations applied to ^He 
p7| . The comparison between our various approaches is interesting, as is that between the 
present data and "exact" results for ^He in Id ||^. Section |VT| summarizes our results and 
conclusions. 



II. HARTREE MODEL CALCULATIONS FOR HELIUM AND HYDROGEN 

Our goal is to determine the threshold coverage for formation of the axial phase. The 
general approach described in this section is analogous to that used in previous studies of the 
problem of layer promotion |]2^j28|| . The assumption is made that, above threshold, atoms 
(or molecules in the case of H2) form two coexisting phases, the axial phase and the shell 
phase. These are characterized by Id densities Na/L and Ng/L, respectively. The condition 
for equilibrium is that the chemical potentials of the two phases coincide: 

f^axial shell (l) 

The axial phase is thought of as a Id phase affected by the "external" Hartree potential 
provided by the shell phase and the host nanotube. Its chemical potential is assumed to 
satisfy 

{Na/ L) — ea + fild{Na/ L) (2) 

The Id approximation ought to be valid here in that we focus on the regime when the rms 
displacement transverse to the axis is small compared to the axial phase interparticle spacing 
and the axial and shell atoms are well separated spatially. Here fiid{Na/L) is the chemical 
potential of a Id fluid at density N^/L and is the eigenvalue of the atoms in the Hartree 
potential. At the threshold for forming the axial phase, fiid assumes its lowest value, the 
ground state cohesive energy of the Id fluid. This energy has been found previously for ^He 
to be of order 2 mK, which is negligible small in the present context, so we ignore it for ^He 



and the other cases studied here 29 



The energy ea was computed from a numerical solution of the Schrodinger equation for 
ground state (zero angular momentum) atomic motion in the Hartree potential: 

Vtot{r) = Vcir) + eVH{r) (3) 

where Vc{r) is the potential near the axis due to the carbon atoms alone computed in the 
assumption of smooth tube walls as in previous studies 0]. The second term is a Hartree 
interaction, where 6 is the density of gas atoms/molecules in the shell, and 6 Vff(r) is the 
interaction due to the shell: 

VH{r) = I rfr'|vl/,(r')pf/,,(|r-r'|) (4) 

where is the wave function of gas atoms/molecules in the shell and Ugg is the gas-gas 

interaction. We make the further approximation of assuming that the atoms in the shell 
phase are narrowly confined in the radial direction: 

|vl/,(r)p ~5(r-i?) (5) 

where R is the radius of the cylindrical surface near which the atoms are situated. Then the 
evaluation of Vh reduces to integrations over the longitudinal and azimuthal coordinates, as 
shown in Ref. |0. 

The ground state eigenfunction of the Schrodinger equation for an atom close to the axis 
is of the form \l/a(r) = f{r) exp{ikz), where r and z are the cylindrical coordinates of the 
atom's position vector r. The eigenvalue is then determined from the differential equation 
satisfied by the ground state radial wave function, /(r): 

We solve this equation numerically. 

Figures |1| and ^ display for the He isotopes the two chemical potentials appearing in 
Eq. ^ their curves' crossing occurs at the threshold at which the axial phase appears. The 
abscissa in these figures is a 2d density of the shell phase: 

9 = Ns/{27tRL) (7) 
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As seen in figures |l| and the eigenvalue ea is a nearly linear function of the shell phase 
density in the neighborhood of the transition, as expected from a perturbation theory of the 
Hartree interaction, based on Eq. |^. 

We have made the simplest plausible assumption for the function UsheiiiQ) - that it 
coincides with the 2d chemical potential on graphite, apart from an additive constant arising 
from the stronger substrate attraction and resulting larger binding energy in the tubes: 

fJ'graphite (9) + AEb (8) 



Here AEf, is the difference in binding energy between the result for the nanotubes [T^ and 
that for graphite |]26| , |30| . The justification for this implicit neglect of the effect of curvature 
appears in the Appendix. Our values of these functions fisheii{9) are taken from theoretical 



studies of Ref. p^ , |25[| which are consistent with experiments reported in Ref. ||26|| . 

The key properties of the shell and axial phases at threshold are displayed in Tables | and 
m and Fig. ^ for the isotopes ^He and ^He, respectively. There are no significant differences 
between the qualitative behaviors of the two isotopes. Consider, for example, the case of 
^He within nanotube having a radius of Rc = 6 A. Note that the He shell radius, R = 3.06 
A, is about 3 A closer to the axis than is the carbon shell; R increases by ca. 1 A as Rc 
increases by 1 A. As Rc increases, the energy levels of the shell (Eq) and axial states at 
threshold increase monotonically. The threshold density 6c of the shell state is insensitive 
to this change, remaining close to the known value 0.115 A^^ for monolayer completion 
density on graphite [^,0. However, the threshold chemical potential fic varies much more; 
it increases by about 20 K for each 1 A change in R; for Rc = 8 A, /ic approaches the value 



-35 K for monolayer completion on graphite . Note that the potential energy responsible 



for the axial state changes by an even larger amount than /ic. The minimum value of the 
total potential energy, Vtot{rmin), changes by 90 K as Rc increases by 2 A, while fic changes 
by only 40 K. The reason is that the zero point kinetic energy also changes significantly, 
as can be understood from the uncertainty principle. The kinetic energy is much larger at 
Rc = 6 A, for which state the root mean square radial coordinate is r^ms = 0.41 A, than at 
Rc = 8 A, for which r^ms = 2 A. Note from the table that the carbon shell's contribution 
to the potential energy Vc{r) is typically 40% of the total potential energy while the He 
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shell contribution is about 60%; these proportions change slowly as Rc changes. Finally, 
we note the qualitative change in the axial state wave function seen in Fig. ^. For the case 
of Rc = 8 A, the axial state's probability density is no longer confined to the immediate 
vicinity of the axis, exhibiting a maximum near r = 2 A. One observable property which is 
determined by this spread is the momentum distribution function, which ought to be much 
narrower in the case of such a dispersed wave function. We intend to consider this topic in 
future work. 

Figures § and |] and Table |Tl| present the results of analogous calculations for H2 adsorp- 
tion in nanotubes of varying size. The data corresponding to Rc = 8 A look qualitatively 
similar to those for He except with different numerical values, of course. As for ^He, the 
threshold values of coverage and chemical potential are similar to those for monolayer com- 
pletion on graphite {6^ = 0.094 A~^, /x^ = -244 K) Q. 

For Rc = 7 A, the axial phase threshold value of the chemical potential for H2, i.e. the 
crossover seen in Fig. |^, moves to a much lower value [fic = —385 K) than for Rc = 8 
A (-261 K). For the smallest case studied, Rc = 6 A, we find no crossing of the chemical 
potential curves. This means simply that nanotubes which are so small do not produce an 
axial phase, according to our calculations. 



III. PATH INTEGRAL SIMULATIONS FOR HYDROGEN 



As a check on the validity of the Hartree calculations we have performed molecular 



simulations of hydrogen adsorption in nanotubes, using the path integral formalism p4 
implemented in the grand canonical ensemble. The algorithm we have used follows the 
previous work of Wang and Johnson ||35|j36[] . The Silvera-Goldman isotropic interaction 
potential was used for the H2-H2 potential. The nanotube-hydrogen potential used in the 
path integral simulations was the same as that used in the Hartree calculations. The path 
integral calculations were carried out at finite temperature, so we do not expect complete 
agreement with the Hartree method. In particular, the shell-axial first order phase change 
will disappear in favor of a continuous transition as the axial phase is populated. Calculations 
were performed at 5 and 10 K for nanotubes of radius 5, 6, 7, and 8 A. In the path integral 
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formalism each quantum molecule is replaced by a classical ring polymer containing some 
number of beads. The accuracy of the simulations depends on the number of beads used 
in each ring, with the results becoming exact as the number of beads per ring approaches 
infinity. We found that 50 beads per ring were sufficient to obtain convergent results at both 
5 and 10 K. 

The shell and axial phase densities of hydrogen in the 7 A tube are plotted in figure |^ 
as a function of the chemical potential. The shell phase density remains fairly constant 
while the axial phase increases dramatically at a chemical potential of about —375 K, in 



reasonable agreement with the Hartree value of /ic = —385 K given in Table |T|. The results 
at a temperature of 5 K (not shown) are very similar, but the population of the axial phase 
commences at a somewhat higher chemical potential. The shell and axial phase densities for 
the 6 A radius nanotube are shown in figure ^ The path integral calculations demonstrate 
that population of the axial phase does occur, with the onset at a chemical potential of 
about —325 K. This is inconsistent with the Hartree calculations, which do not show the 
formation of an axial phase at any chemical potential for this small value of R. 

Figure ^ is a plot of the density profile in the 6 A nanotube at three different chemical 
potentials. Note from figure ^ that the location of the shell phase is pushed progressively 
outward in response to the increased population of the axial phase. The disagreement 
between the Hartree and path integral calculations can be ascribed to the fact that the 
Hartree model does not account for the shift in the location of the shell phase due to the 
presence of the axial phase, and therefore the Hartree model misses the axial phase observed 
in the simulations. The density profile for the 7 A tube is plotted in figure ^ for a series of 
chemical potentials. Note that that shell phase is essentially unperturbed by the presence 
of the axial phase in the 7 A nanotube. 

The coverage and density profiles for the 8 A nanotube are plotted in figures ^ and 



11, respectively. Note that the onset of axial phase population occurs at a higher chemical 



potential (larger bulk pressure) than for either the 6 or 7 A nanotube. This is somewhat 
surprising given that the shell phase must be compressed in order to make way for the axial 
phase in the 6 A tube. The onset of population of the axial phase in the 8 A nanotube 
occurs at about n = —298 K, in fair agreement with the Hartree value of Hc = —261 K from 
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Table |T|. Note from figure [TT| that it is evident that the 8 A nanotube does not actually 
have an axial phase, but rather a second shell phase, centered about 2 A from the center of 
the tube. This agrees with the Hartree calculations of figure |[ 



IV. HELIUM MIXTURES 

We now consider the case of isotopic mixtures of helium adsorbed at medium to high 
density within the tube. In the analogous situation on the graphite basal plane, there is 
complete isotopic separation at T = We expect the same behavior in the nanotubes 



at coverages below that required for producing the axial phase. In this case, the criterion 
determining the distribution of the two phases on the surface is equality of their spreading 



pressures. Figure 12 presents the coexisting 2d densities of these phases, as derived from 



application of this criterion to the equation of state data of Bruch |^9| and Figure |T3| shows 
a schematic view of the mixture configuration, determined by computing the threshold for 
the axial phase. We found that ^He atoms first go to the axial region surrounded by the ^He 
atoms at virtually the same spreading pressure at which ^He atoms begin to go on the axis. 
The preferential binding in the ^He region is a consequence of the higher 2d shell density 
at ^He-'^He coexistence. Our numerical results for this onset condition appear in Figure H 



and Table IV 



V. DENSITY FUNCTIONAL STUDY OF ^HE 

We consider that the energy of an ensemble of helium atoms in the nanotube has the 
form 

E[p]=Eo[p]+ J dvp{r) Vc (r) (9) 

where Vc{r) is the adsorption potential due to all the carbon atoms in the nanotube evaluated 
at point r and Eq[p] is the Finite Range Density Functional (FRDF) that represents the 
energy of inhomogeneous ^He. The form of Eq[p] and the discussion of formalism can be 



found in Ref. [B^. Such density functional have previously been applied to the study of 



adsorbed films |^0|] and impurity atoms in liquid He. |^ 
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At zero temperature, the density of the hehum atoms is p(r) = A^|\l/(r)p, where 
\&(r) is the single particle ground state wave function that minimizes the energy E[p\. The 
minimization procedure leads to the following Hartree-Fock equation for 

(-(nV2m) V^ + t/ [p]) ^ = /i ^ (10) 

Here U[p\ is the effective single particle interaction and is equal to 

U[p] = SE/5p. (11) 

The ground state wave function, \E', depends only on the distance to the axis of the nanotube, 
r, and equation ^ becomes one- dimensional. We compute \E' and /i by solving equation |0 
self-consistently. The results are displayed in figures p!5|-|l7| , for the 6 A nanotube. 



In figure |T5| we plot the density as a function of the radial coordinate r for different 
values of the chemical potential. For small p the density has only one peak centered at 
R = 3.06 A. As long as the value of p increases, the height of this peak grows, reaching 
densities much larger than the saturation density of bulk '^He (0.022 A^'^), suggesting the 
formation of a solid phase (which is not accurately described by our method). For larger 
fi the axial peak becomes appreciable and continues to grow as fi increases, while the shell 
peak remains constant. The picture agrees very well with the model we proposed in section 
1 

We evaluated the number of atoms in the axial region, Na, and in the shell region, A^^, 



by integrating each of the peaks of the density. In figure |Tg we plot the axial density Na/L 
and the shell density 6 = Ns/ {L2ttR) vs. /i. As we can see Na/L is non-zero for p greater 
than pc = -108 K. This axial phase threshold corresponds to N/ L=2.2 A^^ and to 0c=O.113 
in agreement with our previous Hartree model's prediction (see table |). 
The energy per atom and the chemical potential /i are displayed in figure 0. Note a kink 
in the /i curve at the axial phase threshold. This effect can be attributed to the presence of 
a "gas-like" axial phase, since dfi/dN is proportional to the inverse of the compressibility. 
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VI. CONCLUSIONS 



In this work we have presented various calculations pertaining to a hypothetical quasi- 
one-dimensional axial phase of quantum fluids adsorbed in nanotubes. Such a phase ought 
to be thermodynamically distinct from the high density cylindrical shell phase at lower 
density. Our calculations use an approximate Hartree model for H2 and the He isotopes. 
Low temperature path integral calculations for H2 are qualitatively consistent with the 
Hartree results, apart from the case of very small R. For ^He, a density functional model 
yields results which are also consistent with those of the Hartree method. We find such an 
axial phase in all of these cases, except for very small tubes containing H2. We note that 
similar behavior is expected for classical fiuids consisting of small atoms, e.g. Ne and Ar. 
^ Both our methodology and the results presented here are qualitatively similar to those 
obtained previously for the problem of determining the monolayer completion of quantum 
gases adsorbed on graphite p4| , |28| . Quantitatively, however, there is a big difference: the 
threshold chemical potential values are lower for small tubes than for a flat surface, due to 
the higher effective coordination of the particles in small tubes. 

We have not addressed here one of the most important questions: what are the properties 
of this axial phase? We expect these to be entirely different from those of the second layer on 
graphite because of the difference in dimensionality. In particular, we anticipate the behavior 
to be that of the so-called "Luttinger liquid" in one dimension. Among the most remarkable 
of these properties for fermions are the absence of a fermi surface and the dominant role of 
spin and density fluctuations at low T. We hope that relevant calculations and experiments 
are soon forthcoming. 
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APPENDIX Curvature correction to energy 



We here compute the leading curvature correction to the potential energy P{R) per atom 
located on a cylindrical surface of infinite length and radius R. It is typical in such cases to 
find an expansion of the form: 



P{R) ^ [ b 



2 



1+^ (12) 



P(oo) V^. 

where 6 is a characteristic length of the system, assumed to be small compared to R in the 
present asymptotic limit. We shall confirm this expectation here. 

It is reasonable to identify the curvature correction as arising primarily from the attrac- 
tive part of the interatomic interaction since that experiences a larger distance scale than 
the repulsive part of the interaction and so becomes comparable to R first as the graphene 
sheet is curved. The attraction varies as 

V{r) ^ (13) 

where C is the interatomic dispersion coefficient. For the case of atoms having 2d density 9, 

-'^^^ = R f dz f d(p , (14) 

ce J J [z^ + {2Rsin{<^/2yf ^ ' 

where the 2 on the left side avoids double-counting and we bear in mind that eventually we 
will need a small distance cutoff (r = a) in order to get a finite answer. The planar limit, 
R = oo, involves only small (p contributions, leading to 

The difference AP between this value and that of P at finite R may be evaluated with an 
asymptotic expansion. We find 



so that the ratio 

AP 3 /a\2 



— 



(17) 



P 16 \RJ 

confirming the anticipated curvature dependence. A plausible near-neighbor cutoff is a = 

(tt/^)^/^ ~ 2 a. For a nanotube of Rc = 7 A, the film is situated at i? ~ 4 A. Then the 
curvature correction to the energy per atom is about 5%. 
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TABLES 

TABLE L Properties of ^He atoms in the axial and shell states at threshold for formation of 
the axial phase, as a function of radius Rc of the nanotube. The quantities r^ms) Vtot, f-min, 
fj-c, R, and Eq are defined in the text. 



Rc (A) 





7 





r (k) 
' rms \-^J 


n 41 






■Vtotirmin) (K) 


166 


97 


79 


-Vc{rmin) (K) 


60 


40 


33 


1"min (-'^) 


0.00 


1.25 


2.41 




0.113 


0.114 


0.117 


-/"c (K) 


98 


75 


54 


-Eo (K) 


197 


180 


170 


R{k) 


3.06 


4.07 


5.07 


TABLE IL Properties 


of ^He atoms in 


the axial and shell states at threshold, 


as in Table |. 




u 


7 


a 
o 


frms 


0.56 


1.03 


1.96 


■Vtotirmin) (K) 


155 


95 


77 


-Vc{rmin) (K) 


60 


40 


32 


rmin (-'^) 


0.00 


1.38 


2.41 




0.110 


0.110 


0.113 


-^ic (K) 


88 


72 


51 


i?(A) 


3.04 


4.05 


5.02 
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TABLE III. Properties of H2 molecules in the axial and shell states at threshold, as in Table |. 



Rc (A) 


6 


7 


8 


f % \ 

1"rms (-^j 




0.55 


1.38 


-Vtotirmin) (K) 


54.0 


437 


311 


-Vc{rmin) (K) 


254 


141 


106 


' min V-^j 


u.uo 


0.44 


1.67 






0.097 


0.099 


-/"c (K) 




385 


261 


-Eo (K) 


749 


689 


651 


R{k) 


2.69 


3.67 


4.65 


TABLE IV. 


Properties of a ^He-'^He mixture 


in the axial and shell states at threshold, 


as in 


Table |. 








Rc (A) 


6 


7 


8 


/ X \ 


0.43 


1.00 


1.92 


-14of(rmm) (K) 


157 


97 


79 


-Vc(rmm) (K) 


60 


32 


19 




0.00 


1.36 


2.35 




0.110 


0.109 


0.112 




0.114 


0.114 


0.117 


-/"c (K) 


88 


72 


51 
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FIGURES 

FIG. 1. Upwards (downwards) sloping curves are the chemical potentials of the cylindrical 
(axial) phases for ^He as a function of 2d shell density, 64, in the cases of nanotubes with radii Rc 
= 6, 7, 8 A, from bottom to top. 

FIG. 2. Same as Fig. ^ but for ^He. 

FIG. 3. Hartree model calculations of probability density (unnormalized) for ^He atoms in 
tubes of radius of 6 A ( — ), 7 A ( ) and 8 A ( — • — ). 

FIG. 4. Curves are as in Fig. |l], except that no curve exists for Rc = 6 A, as discussed in text. 

FIG. 5. Hartree model calculations, discussed in text, of probability density for II2 in tubes of 

radius Rc'= 6 A ( — ), 7 A ( ) and 8 A ( — ■ — ). Note the absence of the axial phase for the 

case of a 6 A tube. 

FIG. 6. Shell density (9) and axial density (Na/L) in a 7 A radius nanotube at a temperature 
of 10 K as computed by path integral grand canonical Monte Carlo. 

FIG. 7. Shell density (9) and axial density {Na/L) in a 6 A radius nanotube at a temperature 
of 10 K as computed by path integral grand canonical Monte Carlo. 

FIG. 8. Density profile in reduced units for hydrogen in a 6 A radius nanotube at 10 K. Note 
that the shell phase is progressively pushed outward as the axial phase is populated. 

FIG. 9. Density profile in reduced units for hydrogen in a 7 A radius nanotube at 10 K. 

FIG. 10. Shell density {9) and axial density {Na/L) in a 8 A radius nanotube at a temperature 
of 10 K as computed by path integral grand canonical Monte Carlo. 

FIG. 11. Density profile in reduced units for hydrogen in a 8 A radius nanotube at 10 K. 
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FIG. 12. Equilibrium densities of the coexisting ^He and ^He phases. 



FIG. 13. Schematic picture showing the separation of He isotopes and ^He atoms (o) moving 
in the ^He (•) regions. 

FIG. 14. The chemical potential of ^He, /X3, as a function of coverage, 64, in the case when the 
axial ^He atoms go within the ^He shell. 

FIG. 15. Density profile of ^He at T = in the case of a nanotube of radius 6 A from FRDF 
calculation. 

FIG. 16. Axial phase total number density and shell density for ^He in the 6 A nanotube, from 
FRDF calculation. 

FIG. 17. Energy per particle and chemical potential for ^He in the 6A nanotube, from FRDF 
calculation. The dashed line shows the value of N/L at which the axial phase appears. 
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